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EXPERIMENTAL TECHNIQUES FOR DETERMINATION OF THE ROLE OF 
DIFFUSION AND CONVECTION IN CRYSTAL GROWTH FROM SOLUTION 

L. Zefiro 

Genoa University Mineralogy Institute 

The first significant studies of the role of the solution around /663 * 
a growing crystal are probably attributable to Lehmann (1877) , who de- 
veloped various techniques for obtaining a qualitative indication of 
the concentration field around a crystal. For example, he utilized 
substances which produced vividly colored colutions, and observed the 
variations in chromatic intensity near the crystal. 

In 1904, Nerst hypothesized that the concentration of the solution 
in contact with the crystal was the saturation concentration, C g . But 
measurement of the index of refraction based on the angle of total in- 
ternal reflection of NaClO^, NaNO^, KA1(S0 4 )2 and 12^0 (as accomplish- 
ed almost simultaneously by Miers) , resulted in a C Q value for the con- 
centration , such that 

C c C n < C 

where is the concentration of the supersaturated solution far away 
from the crystal. 


Subsequently, Berthoud (1912) , Valeton (1924) and Friedel (1926) 
distinguished two phases in the growth mechanism: 

a) material input by diffusion, and 

b) incorporation into the crystal, 

governed by the following equations: /664 

dm/dt = (D/d) {C a - C 0 ) (1) 

dm/dt = K(C Q - C g ) (2) 

where dm/dt = the quantity of solute at the moment when it is being 

augmented by diffusion and is being diminished by crys- 
tallization; 


♦Numbers in the margin indicate pagination in the Italian text. 
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/•* 


D = the diffusion coefficient; 

d = the width of the layer of solution involved in the diffu- 
sion process; 

k = the growth rate coefficient of one face of the crystal. 

By eliminating Cg, which appears in both equations (1) and (2), 
we obtain the following equation: 

d ” /dt * 4 - D (b/K T (C - - V < 3 > 

where — 3 — ~Tn 7 v\ ' = the effective growth rate coefficient, according 
d - (D/K) fco Friedel# 


Berg (1938) , Bunn (1949) and Humphreys -Owen (1949) measured the 
concentration by means of an interferometric technique consisting of 
placing a previously formed NaClO^ crystal in solution between two 
glass plates which have been silvered so as to be partially reflective, 
and which form a small angle (see Figure 1) . Upon illumination by /665 

monochromatic radiation, a system of interference fringes is formed. 

If the planarity of the plates is accurate, and if the index of refrac- 
tion of the solution within the angle is uniform, then these fringes 
are rectilinear. 
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Figure 1. Diagrammatic representation of the apparatus 

utilized by Humphreys-Owen ( Proc.Roy.Soc . , 1949, 
page 220 ) . ~ 
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Local variations in the index of refraction, due to diminished 
concentration near the growing crystal, cause the fringes to curve 
(see Figure 2) . This phenomenon makes it possible to evaluate the 
distribution of the index of refraction, and also the concentration 
distribution (see Figure 3) , based on the linear relationship between 
these two factors. 

Based on the data obtained, it was concluded that the concentra- 
tion was not uniform along the faces of the crystal but was minimal 
at the center of the faces, where, on the other hand, the concentra- 
tion gradient normal to the face was maximal, i.e., proportional to 
the crystallization flow under the hypothesis that the solute incor- 
poration nechanism is purely diffusive. 



Figure 2. Interference fringes 
obtained by means of 
the apparatus shown 
in Figure 1. 



Figure 5. Dendritic growth. 


Figure 3. Concentration distribu- 
tion around a growing 
crystal, obtained by ex- 
amination of an interfer- 
ometric pattern such as 
that shown in Figure 2. 



Figure 4. Normal growth. 
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In order to obtain an immediate qualitative interpretation of the 
fringes, both Humphreys -Owen and, later, Kceuger and Miller (1953) ar- 
ranged the semi-silvered plates in a parallel configuration, and obtain- 
ed the same type of concentration distribution. 

However, this method has the disadvantage that the growing crystal 
does not always completely fill up the space between the two plates. 

Thus, diffusion occurring relative to the upper face can adversely af- 
fect the results. 

From the fact that the face grows planimetrically , Berg deduced 
that there should be a transport of material along the face in one lay- 
er, perhaps too slight to be observable by means of the interferometry : 
technique utilized, comprising a transition between the solid and li- 
quid states, similar to that introduced by Vollmer (1932) for melting 
growth . 

During the 1949 Faraday Society conference on crystal growth, /666 

Frank opposed this hypothesis on the ground of considerations concern- 
ing the average free exchange of solute molecules. He voiced the opin- 
ion that, in relation to the geometry of the system, it was not justi- 
fiable to leave the convection factor out of consideration. The pres- 
sence of convection would of course invalidate the assumption that the 
solute flow was directly proportional to the normal component of the 
concentration gradient. 

Study of solution growth was advanced by a group of researchers 
at the University of Strasbourg: Goldsztaub and Kern (1953), Folleni- 

us (1959) , Quivy (1965) , Itti (1966) and Goldsztaub, Itti and Mussard 
(1969) , all of whom utilized polarized white-light interferometric tech- 
niques based on a Baker interferential microscope or a dual-refracting 
Frangon interferometer. 

In this type of apparatus, two beams interfere so as to create a 
birefringent beam which passes through the specimen. Thus it is poss- 
ible to measure the optical delay between the perturbed zone of the 
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solution and a second zone which is far away enough from the crystal to 
be considered "undisturbed." 

Once the specimen is aligned in a parallel configuration, the in- 
terference fringes around the crystal are equiconcentrically curved with 
respect to the solution. Based on the forms of the fringes, it is poss- 
ible to distinguish two types of growth: 

a) regular growth, in which the fringes, which are approximately 
parallel to the faces, intersect near the corners (Figure 4); 

b) dendritic (or skeletal, depending on the nomenclature adopted) 
growth, in which the fringes follow the contour of the crystal, 
but without ever intersecting (Figure 5) . 

The results obtained have been interpreted to a first approximation, 
taking into consideration diffusion alone, and considering the problem 
to be two-dimensional. Thus, Fick's equation is utilized: 

d 2 c b 2 c 1 . be ... 

— o ~ ~ — (4) 

Sx z dy D dt 

In order to integrate Fick's equation (4), the initial condition is pos- 
ited that at time t = 0, the concentration is constant, 

c (x,y , 0) = c 

00 

as are the limit conditions /667 


c(»,t) = C. 

and the conditions along the face, according to which the growth is 
classified as either dendritic or regular, respectively: 

Cq = constant, or <^c/dn - constant 

Full integration, as accomplished by Boscher, indicates qualita- 
tive agreement between the c(x,y,t) function and experimentally deter 
mined fringe development. 
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The next step consists of comparing the growth rate V m measured 
at the center of the face and the growth rate V c calculated according 
to the following relation: 


v = dm = D ,3c. 
c psdt p l 5x ; 


(5) 


where p is the density and S is the area. 


The solution of the one-dimensional Fick equation, 


3 c 1 3c 


3x 


1 = d * 5t 


( 6 ) 


with suitable contour conditions, derives from the following relation: 


c (x , t ) = C Q - (C. - C Q ) erf ( ) 

2/Dt 


(7) 


where 


2 _ u 2 

erf (—2—) = -f- I e u du 
2/Dt /? 


/ 


The concentration gradient calculated for x = 0 is; 


,3c x C « " C 0 

( 55T ) = " 


'TT Dt 


( 8 ) 


Thus, it is possible to calculate V and to note that in many instan- 

G 

ces 


V - 2V_ 
m c 


(9) 


The factor of 2 can be explained by the fact chat the measured 
growth rate is an average rate, whereas the calculated growth rate is 
an instantaneous Mte. Consideration of the mean value of V c during 
the time interval between T^ and T ^ yields the following relation: 
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V_ = 


To - T, 


5 | 

P l 


,3c. 

5x } 


O T 
12 


,3c. 

( 3F 


°1 T 1 



Since it has been found experimentally that, except for during 
the initial period, (3c/3x) Q remains nearly constant, it can be pos- 
ited : 


,3c» 

( 5F 


°l T i 


" { 5F 


° 1 T 2 


( 11 ) 


whereby : 


V 

c 


2 § «K> 


( 12 ) 


It must be noted that the apparent relation between V and V 

cm 

derives solely from application of the relation in equation (11) , 
which cannot be justified on the basis of diffusion alone. Therefore, 
it is necessary to take into account convection as well. 

On the simplified hypothesis that the convection rate is constant 
with respect to the variation of x and t, the equation governing the 
transport of material becomes : 


dc 

5t " 



+ Vi 


,3c 

5x 


(13) 


where V is the convection rate. For the concentration gradient on 
the face of the crystal, one obtains the following relation: 


,3c. 

( 5F 



(14) 


By assigning increasing values to the convection rate, the curves 
illustrated in Figure 6 are obtained. It would appear that the intro- 
duction c * convection causes a decrease in deviation from the experi- 
mental curve. 
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Figure 6. Variation in growth rate as a function of time. 

1 = theoretical flow for diffusion alone 

2 through 6 = introduction of convection into 

growth rate calculations 
7 = experimental curve 


Study of Growth from Solution by Means of Holographic Interferometry 

Holographic interferometry (which had previously been utilized in 
our laboratory in order to obtain a record of the microscopic topogra- 
phy of the surfaces being studied) was applied to the study of growth 
from solution. 

An advantage of holographic interferometry is that it effects a 
correlation between the two wave fronts and not between the front of 
the wave being studied and a plane wave, as does conventional inter- 
ferometry. 

Furthermore, the coherence of the laser source utilized makes it 
possible to study larger volumes of solution, and thus makes it easier 
to follow the convective phenomena. 

An example of the experimental apparatus assembly utilized is 
shown in Figure 7, and illustrated schematically in Figure 8. 
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Figure 7. Experimental apparatus assembly utilized to study 

growth from solution by means of holographic inter- 
ferometry. 



Figure 8. Schematic diagram of the apparatus shown in Figure 7 

BS = beam splitter L = lenses 

P = polarizers M = mirrors 

0 = objectives H = holographic plate 


The two collimated laser beams, obtained by means of the beam 
splitter BS, are enlarged by objectives 0 and lenses L. One of the 
two beams passes through the specimen, and the other serves as the 
reference beam. Both beams are reflected by mirrors M onto the holo- 
jraphic plate H where they interfere, creating a fringe pattern with 
sinusoidal amplitude transmission and a fringe period on the order 
of one micron, which can be recorded accurately by the holographic 
emulsion, due to the latter's high resolution (2500 lines/mm). 

After development, the plate is replaced exactly in its initial 
position and illuminated by the reference beam. The microfringe pat- 
tern registered on the plate reveals (by means of diffraction) a wave 
front equal to that which has passed through the specimen. 

If there has been any change in the specimen, the wave front 
which passes through it will be affected, and will interfere with the 
wave front reconstructed by the hologram, thus creating interference 
fringes which indicate where an optical change in the specimen has oc- 
curred. 

This technique, when utilized to follow the evolution of the ph 
nomenon, is termed "real-time" interferometry. 

The image of the specimen is filmed by a television camera, record 
ed by a videorecorder and displayed on a monitor. 



Figure 9. Diffusion fringes in a KC1 solution. 
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The first application of this technique is the determination of 
the diffusion coefficient for a layer of solution whose concentration 
is Cq at a height at whi.h it diffuses toward an overlying solvent col 
umn whose height is b - a. 


In order to obtain the concentration distribution c(x,t) at level 
x and tinu t, it is necessary to solve the Pick equation for the fol- 
lowing initial conditions: 


(x,0) = Cq for ■> < x < a 
(x,0) = 0 for a < x < b 

which results in a Fourier series of the following type: 


c(x,t) = C, 


a 

E 


w 

IE 

n=l 


sin(^) 


cos (S**) 


-21 2 Dt 
„ b 


n 


i 

(15) 


The variations in the index of refraction corresponding to changes /671 
in concentration are indicated holographically by means of a pattern of 
horizontal fringes (see Figure 9) . 

Based on the correspondence between c(x,t )/C q and the expression /672 
[n(x,t) - nf^on/n^ where n (| is the index of refraction of the solu- 
tion having concentration C n , it is possible to obtain a reasonable 

-5 2 U 

value for h of 2 x 10 cm /sec. for diffusion of a 10% KC1 solution. 

The growth of NaClO^ in a vessel 20 mm in diameter and 1 mm wide 
was also studied. The vessel was placed in region "A" of the diagram 
shown in Figure 8, in horizontal and vertical positions, the better 
to demonstrate convective movements, with the results shovm in Figure 
10 and Figure 11, respectively. 

Examination of Figure 10 will reveal that the concentration is con- 
stant parallel to each face, except near the corner •- Simply by count- 

ing the number of fringes and referring to the concentration of the 
supersaturated solution at a distance from the crystal, it can be deduc- 
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Figure 10. Growth of an NaClO-j crystal: horizontal view of 

the vessel. 
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Figure 11. Growth of an NaCIO, crystal: vertical view of 

the vessel. 






ed that Cq, the concentration of the solution in contact with the crys- 
tal , assumed a different value depending on the face being studied , and 
has a maximum value in relation to the corner. 

In Figure 11, it will be noted that, due to the established convec- 
tive regime, an ascending "chimney" forms in relation to the upper hori- 
zontal face of the crystal. This "chimney" allows evacuation of the less 
supersaturated solution from the zone nearest the crystal. The other, 
more supersaturated solution replaces it, thus increasing the growth rate. 

3 

Tests carried out on solution volumes on the order of 30 cm , suit- 
ably thermostat-controlled with regard to temperature, have also indicat- 
ed the predominant role played by convection. 

In order to obtain a quantitative interpretation of the experimental 
results, it will be necessary to work in two different directions: 

1) to study capillary gio* th as Marsiglia B. Simone of the Crystal 
Growth Mechanism Research Center is doing. If the crystal occu- 
pies the entire section of the vertically positioned capillary, 
then it would appear to be correct to state that the growth of 
the lower horizontal face is governed by a purely diffusive mech- 
anism. The advantages of the analytic interpretation are oppos- 
ed by the geometry of the problem, which implies an indirect re- 
lationship between the fringes and the concentrations. 

2) to replace the television camera with a quadratic photodiode ma- 
trix, digitalize the data obtained and analyze them by means of 
a computer. 

It can be hoped that in this way it will be possible to describe the 
c(x,y,t) function by taking into account both the diffusion and convection 
mechanisms. 
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